Journal of Experimental Medicine
● Rockefeller University Press
Preprints posted in the last 7 days, ranked by how well they match Journal of Experimental Medicine's content profile, based on 119 papers previously published here. The average preprint has a 0.11% match score for this journal, so anything above that is already an above-average fit.
Ding, H.; Cox, M.; Qiu, Q.; Mirabadi, N.; Bobrovnikov, D.; Gerbino, K.; Scherler, K.; Walters, K.; O'Keefe, G.; Brakenridge, S.
Show abstract
BACKGROUND Traditional organ dysfunction classifications after severe trauma rely on binary categories (e.g., chronic critical illness [CCI] versus rapid recovery), which likely oversimplify the true heterogeneity of postinjury organ dysfunction trajectories and their underlying immune-biologic drivers. We hypothesized that severely injured trauma patients demonstrate multiple, distinct organ dysfunction phenotypes distinguished by unique immune signatures. METHODS We performed a prospective, longitudinal, observational cohort study of severely injured trauma patients at high risk for multiple organ dysfunction at a single Level 1 trauma center. Organ dysfunction was quantified daily through 28 days, death, or discharge using a trauma-specific modified Sequential Organ Failure Assessment (mSOFA) score excluding the neurologic subdomain to isolate systemic organ dysfunction from traumatic brain injury. We utilized group-based trajectory modeling to identify longitudinal mSOFA phenotypes. Plasma immune protein trajectories were compared across phenotypes using mixed-effects models with time interaction. RESULTS Sixty-nine patients were enrolled. Trajectory modeling identified five mSOFA phenotypes (c1-c5) with differences in resuscitation requirements, shock severity, and 28-day mortality, but not baseline demographics, injury severity, or early physiologic parameters. The phenotypes revealed organ dysfunction heterogeneity not fully captured by CCI/rapid recovery classification. Forty of 92 measured immune markers showed phenotype-differentiated trajectories (FDR<0.1). The most fulminant phenotype (c5) was distinguished by IL-17C, CXCL5, and CXCL6, implicating epithelial-barrier and neutrophil-trafficking mechanisms. The phenotype with highest CCI incidence (c3) was distinguished by rising CXCL9, ARTN, and TWEAK, implicating Th1 polarization and tissue-inflammation/repair signaling. The clinically distinct, older-age phenotype (c4) carried the largest differentiating signature (29 proteins), spanning renal injury/repair, endothelial activation, and adaptive lymphocyte modulation. CONCLUSIONS Organ dysfunction after severe trauma is more biologically heterogeneous than traditional clinical classifications suggest. Distinct mSOFA trajectory phenotypes correspond to distinct immune signatures, supporting separable pathobiology rather than a single axis of injury severity, with implications for phenotype-targeted precision immunomodulatory approaches in critical illness.
Kulshrestha, A.; Garrett, M. E.; Nouraie, S. M.; Machado, R. F.; Gordeuk, V. R.; Gladwin, M. T.; Telen, M. J.; Ashley-Koch, A. E.; NHLBI Trans-Omics for Precision Medicine (TOPMed),
Show abstract
Sickle cell disease (SCD) results from homozygosity or compound heterozygosity for a single point mutation in the beta hemoglobin gene. Despite carrying the same underlying causal mutation, patients exhibit highly heterogeneous clinical outcomes. Acute chest syndrome (ACS) is one of the leading causes of mortality in SCD. In this study, we aimed to identify genetic variants associated with ACS occurrence and recurrence (>5 episodes) in two SCD cohorts (OMG-SCD and Walk-PHaSST) of predominantly African ancestry and characterize their downstream molecular consequences. Although none of the SNPs reached the genome-wide threshold in the occurrence meta-analysis, we identified 40 independent genetic loci with p <10-5. However, the recurrence GWAS in OMG-SCD, which compared individuals with no prior history of ACS to individuals with a history of >5 episodes, identified one genome-wide significant locus along with 38 other independent loci with p <10-5. Survival analysis suggested lower survival probability among young ACS recurrence cases compared to controls or occurrence cases. Furthermore, molecular QTL analysis of ACS recurrence identified 42 SNP-protein and eight SNP-metabolite associations with p <10-3. Pathway and enrichment analysis identified secretory, platelet/hemostatic, and complement-mediated pathways, confirming the complex and heterogeneous pathophysiology of ACS. This study highlights the importance of a well-defined phenotype for identification of risk loci and elucidates the complex interplay of hemolytic, inflammatory, and cellular organization pathways that warrants further evaluation through large-scale multi-omic studies of SCD patients. Ultimately, the more extreme phenotype of ACS recurrence may be a more informative phenotype for analysis, risk prediction, and therapeutic targeting.
Kuramoto, K.; Ullah, F.; Liang, P.; Ramirez, I. M.; AlAbdi, L.; Feng, B.; Sakurai, M.; Kang, M.-S.; Dong, H.; Oyeyemi, M.; Qin, Z.; Coghlan, M. A.; Sharma, A.; Mercier, A.; McRae, A.; Belmadani, A.; Contractor, A.; Abuyousef, O.; Hashem, M. O.; Alhashem, A.; Maddirevula, S.; Alkuraya, F.; He, Y.; Davis, E. E.; He, C.
Show abstract
Genetic regulation of neurodevelopmental disorders that impair brain development and social interaction is complex and poorly understood. Although emerging evidence suggests a link between neurodevelopmental disorders and autophagy deficits, congenital autophagy gene disorders are rare in humans, and underlying molecular mechanisms remain unknown. Here we demonstrate that the autophagy protein BECN1 is a regulator of neuronal and social development in vertebrate species. We identify rare deleterious biallelic BECN1 variants that cause neurodevelopmental disorders in two unrelated families. Becn1 disruption also causes neuronal and social developmental defects in mouse and zebrafish models. BECN1 serves as an adaptor to recruit Arc to the E3 ubiquitin ligase UBE3A for polyubiquitination and degradation, thereby maintaining postsynaptic AMPA receptor (AMPAR) stability, AMPAR-mediated excitatory neurotransmission, and normal social interaction. Thus, our study illustrates new biology and a non-canonical mechanism by which BECN1 regulates neurodevelopment in humans, mice and zebrafish, via a crosstalk with the ubiquitin-proteasomal machinery.
Sezin, T.; Oktafina Prayitno, N.; Gupta, Y.; Mousavi, S.; Tang, X.; Bivona, D. L.; Zirpel, H.; Staender, S.; De Luca, D. A.; Papara, C.; Pommerien, L.; Thaci, D.
Show abstract
How inflammation is established and resolved within tissues remains a fundamental question in human biology. Psoriasis is a chronic inflammatory skin disease driven by the IL-23/IL-17 pathway, with IL-17 pathway targeting biologics producing more rapid clinical improvement than IL-23 inhibitors. This difference raises the question of how IL-17 and IL-23 inhibition differentially remodel the cellular and transcriptional landscape of psoriatic skin during resolution. Here, we generated a longitudinal spatial atlas of plaque psoriasis using 2 micron high-definition spatial transcriptomics (HD-SPT), profiling 42 matched biopsies from 11 biologic treatment-naive patients before and during IL-17- or IL-23-targeted therapy. By integrating HD-SPT with a single-cell RNA-sequencing (sc-RNAseq) reference atlas, we resolved the cellular composition and spatial organization of psoriatic skin and mapped its remodeling during treatment. We identified a distinct cluster of SPRR2F+ keratinocytes (KCs) that was spatially restricted to the upper stratum granulosum at the transition to the stratum corneum and was strongly associated with inflamed psoriatic epidermis. SPRR2F+ KCs were enriched for antimicrobial and inflammatory barrier genes, including S100A8, S100A9, LCN2, IL36G, and PI3, and their abundance correlated with disease activity during treatment. In psoriasis patient-derived human epidermal KCs (P-HEK), SPRR2F knockdown attenuated inflammatory and antimicrobial transcriptional responses, supporting a functional role for SPRR2F in sustaining the inflammatory KC response in psoriasis. Longitudinal profiling further revealed distinct trajectories of tissue remodeling following IL-17 and IL-23 inhibition. Following IL-17 inhibition, SPRR2F+ KCs and their characteristic transcriptional signature resolved rapidly, whereas IL-23 inhibitors produced more gradual resolution of KC-centered transcriptional changes accompanied by broader immune and stromal remodeling. Together, these findings established the first HD-SPT atlas of a human skin disease and revealed distinct cellular and transcriptional trajectories of psoriatic skin during therapeutic resolution, identifying SPRR2F as a regulator of persistent KC inflammation.
Chen, E.; Gabernet, G.; Corry, D. B.; Kheradmand, F.; Baden, L. R.; Sekaly, R.-P.; McComsey, G. A.; Haddad, E. K.; Cairns, C. B.; Rouphael, N.; Fernandez-Sesma, A.; Simon, V.; Metcalf, J. P.; Agudelo Higuita, N. I.; Hough, C. L.; Messer, W. B.; Davis, M. M.; Nadeau, K. C.; Pulendran, B.; Maecker, H. T.; Kraft, M.; Bime, C.; Reed, E. F.; Schaenman, J.; Erle, D. J.; Calfee, C. S.; Atkinson, M. A.; Brakenridge, S. C.; Melamed, E.; Ehrlich, L. I. R.; Shaw, A. C.; Hafler, D. A.; Augustine, A. D.; Becker, P. M.; Peters, B.; Ozonoff, A.; Diray-Arce, J.; Levy, O.; Bosinger, S. E.; Eckalbar, W.; IMPACC
Show abstract
Host genetic variation influences heterogeneity in infectious disease severity. Genome-wide association studies (GWAS) have identified 69 lead single nucleotide polymorphisms (SNPs) associated with critical COVID-19, but most lack known regulatory mechanisms. To date, the genetic cis-regulatory effects of critical COVID-19 GWAS lead SNPs have not been systematically characterized across acute illness and convalescent phases in a COVID-19 specific cohort. Here, performing severity-stratified, longitudinal cis-expression quantitative trait loci (cis-eQTL) mapping and colocalization with critical COVID-19 GWAS signals, we find that host cis-regulatory effects vary across clinical severity and disease phase. We analyzed peripheral blood mononuclear cells from 1,031 patients in the Immunophenotyping Assessment in a COVID-19 Cohort (IMPACC), a prospective multicenter study of adults hospitalized with COVID-19, profiled across 10 timepoints spanning acute illness and convalescence. We identified 23 GWAS lead SNPs that colocalized with cis-eQTLs at 26 candidate causal genes, with 16 detectable only when severity stratification and longitudinal sampling were incorporated. A colocalized, protective allele (chr19:4,717,660) at the TMIGD2 locus upregulated TMIGD2 expression exclusively during the acute phase in the prolonged-critical/fatal cohort, with no detectable effect in the rapid-recovery cohort. Our results show that host genetic effects on critical COVID-19 can be context-dependent, with cis-regulatory effects varying across clinical severity and disease phase. Severity-stratified, longitudinal cis-eQTL designs have greater sensitivity to identify candidate causal genes for critical COVID-19 GWAS associations than single-timepoint cross-sectional studies.
Torreni, L.; Schluter, A.; Marti, S.; Ruiz, M.; Velez Santamaria, V.; Rodriguez-Palmero, A.; Artuch, A.; Guilera, C.; Arias, A.; Mochel, F.; Aguilera-Albesa, S.; Cantarin, V.; Casasnovas, C.; Artuch, R.; Koehler, W.; Pujol, A.; FOURCADE, S.
Show abstract
X-linked adrenoleukodystrophy (X-ALD) is the most common monogenic disorder affecting brain white matter, presenting a broad clinical spectrum from slowly progressive adrenomyeloneuropathy (AMN) to rapidly advancing cerebral demyelination. Factors influencing the variable disease expressivity remain poorly understood. We conducted shotgun metagenomic profiling of fecal samples from X-ALD patients and integrated these data with targeted plasma metabolomic analyses of amino acids and key metabolic intermediates. We identified distinct microbiome-associated metabolic signatures across X-ALD phenotypes. Notably, Bacteroides cellulosilyticus and Helicobacter apodemus abundance were significantly altered in AMN compared with childhood CALD patients. In parallel, X-ALD patients exhibited evidence of microbial metabolic remodeling characterized by altered amino acid metabolism, reduced folate biosynthetic potential, and enrichment of stress-adaptive pathways, including redox-associated functions, glyoxylate shunt activity, lipid A modification, lipopolysaccharide assembly, biofilm adhesin biosynthesis, and curli fiber production. Consistent with these microbial changes, plasma levels of branched-chain amino acids and methionine, as well as serum folate concentrations, were significantly reduced. Importantly, methionine levels showed an inverse correlation with neurofilament light chain (NfL) concentrations in patients with cerebral disease. Our findings reveal a previously unrecognized link between gut microbial metabolic capacity and systemic metabolic alterations in X-ALD. Our results reinforce previous evidence implicating methionine metabolism in disease severity. Collectively, these data support a modulatory role of the microbiota-gut-brain axis in neuroinflammatory and mitochondrial-related processes underlying X-ALD progression and highlight new opportunities for biomarker development and microbiome-informed therapeutic interventions.
Chaves Filho, A. M.; Nascimento, E. F.; Goncalves, L. D. C.; Dias Filho, F. J.; Oliveira, T. Q.; Macedo, D. S.; Dal-Pizzol, F.; de Moraes, M. E. A.; Barichello, T.
Show abstract
OBJECTIVE: Sepsis-associated encephalopathy (SAE) and sepsis-associated delirium (SAD) lack confirmatory laboratory tests or biomarkers. We evaluated whether the Th17/interleukin (IL)-17 axis or neuroglial-injury biomarkers reflect acute brain dysfunction beyond systemic inflammation and whether biomarker profiles differ between these phenotypes. DATA SOURCES: EMBASE, MEDLINE/PubMed, and the Cochrane Library from inception to August 1, 2026, supplemented by backward and forward citation searching, following PRISMA 2020. STUDY SELECTION: Observational studies of adults with sepsis classified by SAE or SAD status and reporting circulating biomarker concentrations by neurological group or biomarker-outcome associations. Eighteen studies (2,388 participants; 1,863 with sepsis) were included. DATA EXTRACTION: Two reviewers independently extracted study, phenotype, assay, and outcome data. Risk of bias was assessed with the Newcastle-Ottawa Scale and certainty with GRADE. DATA SYNTHESIS: Random-effects models pooled standardized mean differences (SMDs) separately against non-septic and septic controls and odds ratios (ORs) for biomarker-outcome associations. IL-6 was the only mediator elevated in both comparative frameworks and associated with outcomes (non-septic controls: SMD, 0.54 [0.25-0.82]; septic controls: SMD, 1.05 [0.51-1.59]; OR, 1.85 [1.50-2.29]). S100B, neuron-specific enolase, IL-10, tumor necrosis factor; and interferon-gamma; were elevated versus septic controls, whereas IL-17 was modestly elevated only versus non-septic controls (SMD, 0.37 [0.01-0.73]). The IL-6 septic-control effect was concentrated in SAE (p<0.0001) and remained robust in leave-one-out and risk-of-bias sensitivity analyses. Certainty was low to very low, and all 14 computable prediction intervals included the null. CONCLUSIONS: IL-6 emerges as the leading circulating candidate biomarker of sepsis-associated brain dysfunction, while S100B and NSE capture neuroglial injury beyond sepsis. Current evidence does not establish a discriminating role for the Th17 axis. Future studies require standardized neurological phenotyping and appropriately matched septic controls to establish diagnostic value.
Taylor, R. W.; Smith, S. J.; Ait-El-Mkadem Saadi, S.; Leighton, J. K.; Taylor, L. S.; Frey, A. M.; Nelson, M. G.; Castaneda-Tamez, P.; Genin, E. C.; Caswell, R.; Dannoura, A.; Shamsheldin, H. E.; Fragaki, K.; Rouzier, C.; He, L.; Hemingbrough, C. V. Y.; Fassad, M. R.; Fouilhoux, A.; Shuaib, T.; Alharby, E.; Rehsi, P.; Rodenburg, R.; Rinne, T.; Janssen, M. C. H.; Hashem, M. O.; Hamid, H.; Bakur, K.; Alrabee, H.; Mohamed, S.; Alansari, R. M.; Al-Owain, M.; Blakely, E. L.; Winter, M.; Almontashiri, N. A. M.; Morris, A. A.; AlAbdi, L.; Alkuraya, F. S.; Robinson, H. K.; McFarland, R.; Alston, C. L.
Show abstract
Mitochondrial disease arises from variants in either the mitochondrial or nuclear genome that impair mitochondrial oxidative phosphorylation. Complex III (ubiquinol:cytochrome c oxidoreductase, CIII2) is a key enzyme of the mitochondrial respiratory chain, catalyzing electron transfer from reduced coenzyme Q10 (ubiquinol) to complex IV, via cytochrome c. CIII2 has one subunit encoded by mitochondrial DNA, the remainder being nuclear encoded. Individuals with CIII2 deficiency present with significant clinical heterogeneity, on account of pathogenic variants in genes encoding structural CIII2 subunits or assembly factors. UQCC1 (GenBank: NM_018244.5) encodes the mitochondrial complex III assembly factor Ubiquinol-cytochrome c reductase complex assembly factor 1, which is currently not associated with human disease. Here we describe nine subjects from eight unrelated families harboring biallelic UQCC1 variants with heterogeneous, early-onset mitochondrial disease presentations and lactic acidosis. Subject-derived fibroblasts demonstrated a decrease in the steady-state levels of UQCC1 and its interacting partner UQCC2, with structural modelling supportive of variant pathogenicity. BN-PAGE and quantitative proteomics demonstrated decreased CIII2 abundance whilst mass spectrometry (MS)-complexome profiling corroborated these findings, identifying impaired respiratory supercomplex formation and a delay in the incorporation of newly synthesized CIII2 subunits into native assemblies. Lentiviral transduction experiments in subject-derived fibroblasts restored the levels of both UQCC1 protein and CIII2 activity. Collectively, these findings provide convergent genetic, biochemical, proteomic and functional evidence that bi-allelic UQCC1 variants disrupt CIII2 assembly, establishing UQCC1 deficiency as a cause of mitochondrial CIII2 deficiency.
Whittle, B. J.; Williams, M.; Lomax-Browne, H.; Sorensen, L.; Maughan, R. T.; Porteous, M.; Chakrabarty, A.; Luqmani, R. A.; Pickering, M. C.; Morgan, A. W.; Peters, J. E.
Show abstract
Giant cell arteritis (GCA) is an immune-mediated vasculitis that can cause serious morbidity including sight loss, stroke and aortic aneurysm. There is an unmet need for blood-based biomarkers and novel therapeutics. We measured plasma protein levels using the SomaScan v4.1 platform, comprising 7,273 protein measurements (SOMAmers) reflecting 6,596 unique protein targets. We profiled the proteome in 113 patients who presented acutely with possible GCA. 56 had biopsy-proven GCA and 57 did not have GCA (non-GCA). We found 144 SOMAmers increased and 151 reduced in GCA. Dysregulated biological pathways included IL6-JAK-STAT3 signalling, chemokines for both myeloid (e.g. CXCL2, CXCL3, CCL18, CCL23) and lymphoid cells (e.g. CXCL9, CXCL13), and the complement system (e.g. C9, CFHR5, C5a, Factor B). Enzyme-Linked Immunosorbent Assays (ELISAs) confirmed elevated plasma levels of C5a, the terminal complement complex, Factor B and CFHR5 in GCA. We developed a 26-protein signature that distinguished GCA from non-GCA (area-under-the curve, AUC, 0.93) and validated the diagnostic signature in external datasets. Our data implicate terminal complement pathway activation in GCA, raising the possibility of complement-directed therapies. We show proof-of-principle that a multiplex protein biomarker could be used as a diagnostic adjunct in suspected GCA.
Robinson, K. R.; Curtis, S. W.; Kurtas, N. E.; Sanchis-Juan, A.; Manchel, A.; Adeyemo, W. L.; Beaty, T. H.; Butali, A.; Buxo Martinez, C. J.; Cutler, D. J.; Epstein, M. P.; Feingold, E.; Gowans, L. J.; Hecht, J. T.; Moreno Uribe, L.; Padilla, C. D.; Shaw, G. M.; Weinberg, S. M.; Murray, J. C.; Cotney, J.; Talkowski, M.; Marazita, M. L.; Brand, H.; Leslie-Clarkson, E. J.
Show abstract
Structural variants (SVs) are larger genetic rearrangements such as deletions, duplications, and insertions and other forms of structural alterations. Given their large size, they can have an outsized effect on phenotype, however, the contribution of SVs to the genetic architecture of orofacial clefts (OFCs) is not well understood. In this study, we sought to better understand the contribution SVs to OFC risk by characterizing de novo SVs (dnSVs) in 1,549 probands with an OFC. Overall, we called 375 dnSVs in 307 probands (0.242 dnSVs per genome). The dnSV rate per genome was highest among individuals with cleft lip (CL, 0.29), followed by cleft lip and palate (CLP, 0.25) and cleft palate only (CP, 0.20), though differences across subtypes were not statistically significant. There were also no differences between sex or syndromic status, but we did identify a significant difference in the rate of duplications across OFC subtypes (p=0.005). Compared to CP (0.028), duplication rates were significantly higher for both CL (0.121; p = 0.003) and CLP (0.096; p = 0.002), suggesting a potential etiologic mechanism for CL/P. We systematically curated each dnSV to classify its potential for influencing OFC risk by examining the genes and regulatory elements disrupted by the SV and evidence for their role in craniofacial development or OFCs. We identified 22 dnSVs in regions associated with recurrent genomic disorders, including seven deletions and one duplication of the well-established OFC risk locus 22q11.2. Additionally, we identified 12 large (>1Mb) SVs, with eight meeting criteria for a likely pathogenic classification based on total gene involvement, and seven dnSVs disrupting well-known genes associated with autosomal dominant forms of OFC. The most notable was SATB2, for which we found four dnSVs, all in individuals with CP. Next, by combining evidence from recurrent dnSVs and de novo single nucleotide variants with Hi-C contact maps and predicted chromatin states from developing craniofacial tissues, we found new candidate genes for OFC risk including FHAD1, NBEA, NR3C2, and SPHKAP. Overall, dnSVs classified as likely pathogenic accounted for putative molecular diagnoses in approximately 1.67% of the cohort, with an additional 0.96% having variants of unknown significance with plausible risk for OFCs. Our results support dnSVs as an important component of OFC etiology and suggest their inclusion in analyses may improve diagnostic yield for affected individuals.
Molotkov, I.; Reeve, M. P.; Kolosov, N.; FinnGen, ; Tsao, H.; Palotie, A.; Ripatti, S.; Daly, M. J.; Gupta, S.; Artomov, M.
Show abstract
Hypothyroidism has been associated with increased skin cancer risk, whereas greater genetic susceptibility to hypothyroidism is associated with lower risk. These opposing associations suggest that the average risk associated with hypothyroidism may conceal substantial heterogeneity among patients, and that genetic background could help identify those at highest risk. In this study, we analyzed more than 650,000 individuals across All of Us and FinnGen cohorts to determine whether a hypothyroidism polygenic risk score (H-PRS) stratifies skin cancer risk among patients with autoimmune hypothyroidism and whether integration of a skin cancer polygenic risk score (SC-PRS) further refines risk. In FinnGen, autoimmune hypothyroidism was associated with increased skin cancer hazard (hazard ratio [HR] = 1.10; 95% CI [1.06, 1.15]), whereas H-PRS was inversely associated with risk (HR=0.97 per 1-SD increase; 95% CI [0.95, 0.98]). The hypothyroidism x H-PRS interaction was also negative (HR=0.94; 95% CI [0.91, 0.98]), indicating greater hypothyroidism-associated risk at lower H-PRS. The hypothyroidism and H-PRS associations replicated in All of Us, and similar H-PRS risk gradients were observed for melanoma, basal cell carcinoma, and squamous cell carcinoma in both cohorts. Among hypothyroid patients, predicted 10-year risk increased with higher SC-PRS and lower H-PRS. Combined clinical-genetic profiles had HRs ranging from 2.32 at the 20th H-PRS/80th SC-PRS percentiles to 3.37 at the 5th/95th percentiles, exceeding published point estimates for the established dermatological risk factors. Among women aged 55, representing the more commonly affected sex at approximately the median age of diagnosis, 10-year risk of skin cancer was 14.9% (95% CI [8.4%, 20.9%]) in the highest-SC-PRS/lowest-H-PRS quartiles versus 2.8% (95% CI [1.0%, 4.7%]) in the opposite quartiles. Overall, these findings support PRS-guided risk stratification for targeted skin cancer prevention and surveillance.
Li, Y.; Zhou, H.; Yeh, F.-Y.; Sun, Z.; Zheng, J.; Alzaidi, S.; Guan, H.; Li, J.; Seger, D. L.; Wang, L.; Li, M.; Li, X.; Chen, T.; Plasek, J. M.; Li, H.; Jiang, A.; Hong, P.; Li, F.; Zhang, W.; Cui, L.; Bates, D. W.; Zhou, L.; He, Y.; Tao, C.
Show abstract
Biomedical artificial intelligence (AI) increasingly relies on heterogeneous data sources, yet differences in data structure, terminology, evidence quality, and provenance complicate their integration into reliable knowledge resources. Here, we develop an evidence-aware framework for processing and integrating heterogeneous pharmacovigilance data from FDA Adverse Event Reporting System (FAERS) reports, biomedical literature, clinical trials, and electronic health records in Mass General Brigham. The framework combines large language model (LLM)-based information extraction, ontology-driven semantic normalization, cross-source evidence assessment, statistical and machine learning-based prioritization, and ingredient-based adverse drug event (ADE) relationship inference to transform heterogeneous evidence into a structured and evidence-based knowledge resource. The constructed knowledge system contains 11,588 drug nodes, 44,843 ingredient nodes, 3,275 ADE nodes, 993 indication nodes, and more than 2.47 million relationships, with each drug-ADE relationship assigned a confidence score. We further developed a drug-level web interface for accessing the knowledge resource and validated its utility in downstream AI tasks. Specifically, incorporating the developed knowledge resource improved the comprehensiveness, relevance, and overall quality of LLM-based drug-specific ADE synthesis and improved the precision and F1 score of clinical ADE relation extraction. These results demonstrate that AI-assisted processing combined with evidence-aware curation can transform heterogeneous biomedical evidence into reusable, AI-ready knowledge with explicit provenance and evidence strength information. The framework provides a general strategy for integrating heterogeneous biomedical data to support reproducible knowledge construction and downstream biomedical AI applications.
Hoffman, K. L.; Weil, A. A.; Wetzler, E. A.; Akinola, B. O.; Beauge, A.; Andy, C.; Adams, A.; Frivold, C.; Cox, S. N.; Ehmen, B.; Hatchie, T.; Clark, E. L.; Drummond, M.; Englund, J.; Han, P. D.; Kuntz, J. L.; Lo, N. C.; Lockwood, C.; Starita, L.; Mularski, R. A.; Naleway, A. L.; Carone, M.; Chu, H. Y.
Show abstract
Background: The relationship between recent viral infection and subsequent susceptibility to, or protection against, other viruses remains unclear, particularly across ages and symptom status. Methods: A prospective cohort study enrolled participants aged 6 months to 49 years in Washington and Oregon, USA, 2022-2024. Weekly nasal swabs were tested for SARS-CoV-2, respiratory syncytial virus (RSV), and influenza A/B; additional swabs were collected during symptomatic illness and were also tested for parainfluenza (PIV) 1-4, endemic human coronaviruses (hCOV), adenovirus, human metapneumovirus, human rhinovirus (hRV), and enteroviruses. Proportional hazards regression models were used to characterize associations between incident symptomatic infection (illness) and history of recent illness due to other viruses (within 1-7, 8-14, 15-30, 31-60, 61-90, 91-180 days), adjusted for geographic region, age, sex, daycare/school/work attendance, household size, and number of children in the household. For SARS-CoV-2, influenza, and RSV, these analyses were repeated to study incidence of infection (including asymptomatic cases). Findings: Of 3,560 participants, 2,969 (83%) had [≥]1 illness over 4,299 person-years. Adjusted hazard ratios (aHRs) in the short term ([≤]2 weeks) after an illness were predominantly below 1, indicating that recent viral illness was often associated with lower risk of subsequent illness with other viruses. Across age groups, recent hRV illness was associated with lower risk of hCOV, PIV, and RSV illness (aHRs 0.07-0.27). SARS-CoV-2 illness was associated with a higher short-term risk of adenovirus and PIV illness (aHRs 3.20-3.33). In children aged 6 months-4 years, recent influenza illness was associated with higher short-term risk of SARS-CoV-2 illness. Including asymptomatic infections generally attenuated these associations. Interpretation: Some viral illnesses may alter short-term susceptibility to subsequent viral illnesses. Interactions may differ by age and symptom status. Funding: US Centers for Disease Control and Prevention; Washington State Department of Health.
Wang, X.; Macieira, T. G. R.; Hong, Y.-R.; Revere, L.; Xie, Z.
Show abstract
Objective To examine the association between patient-reported shared decision-making (SDM) and documented advance care planning (ACP) among adults with Alzheimer disease and related dementias (ADRD), and to assess whether this association differed between adults with and without ADRD. Design Retrospective observational study using linked survey and electronic health record (EHR) data. Setting and Participants Participants aged [≥]50 years in the All of Us Research Program who completed a patient-reported SDM measure and had linked EHR records. The primary analytic cohort included 3,635 adults with ADRD; 123,681 adults without ADRD comprised a secondary comparison cohort. Methods SDM was measured by how often clinicians asked participants about their opinions or preferences regarding medical care or treatment and was categorized as high (always/most of the time) or low (some of the time/ None of the time). Documented ACP was identified within 3 years before or after the SDM assessment. Modified Poisson regression with robust variance estimated prevalence ratios (PRs) and adjusted PRs (aPRs). Results Among adults with ADRD, 60.6% reported high SDM. Documented ACP prevalence was 5.1% among participants with high SDM and 3.6% among those with low SDM (P=0.027). High SDM was associated with greater ACP prevalence in the unadjusted model (PR, 1.44; 95% CI, 1.04-1.99); the association remained positive but was attenuated after full adjustment (aPR, 1.28; 95% CI, 0.93-1.77). Among adults without ADRD, high SDM was associated with documented ACP (aPR, 1.17; 95% CI, 1.03-1.33). Conclusions and Implications Among adults with ADRD, high SDM was associated with greater ACP prevalence in unadjusted analyses, although the fully adjusted association was not statistically significant. In contrast, higher SDM remained significantly associated with ACP among those without ADRD. The low prevalence of ACP highlights an important opportunity to expand discussions of current and future care preferences as a routine component of patient-centered care.
Lombardi, G.; Leung, C. C.-Y.; Burge, R.; Zelaya, F.; Bhattacharyya, S.; Velayudhan, L.
Show abstract
Neuropsychiatric symptoms (NPS) are highly prevalent in Mild Cognitive Impairment (MCI), conceptualized as an intermediate clinical state between normal aging and dementia. Although neuroimaging studies have reported associations between NPS and brain alterations in MCI, longitudinal evidence on NPS and their relationship with multimodal MRI measures remains limited. We aimed to investigate (1) longitudinal changes in NPS, volumetric assessment and measurement of regional cerebral blood flow (CBF) within the medial temporal lobe (MTL); (2) associations between longitudinal changes in MRI measures and NPS; (3) whether baseline MRI measures predicted longitudinal changes in NPS. Thirty individuals with MCI aged over 65 years were recruited, of whom eighteen followed up for a mean period of 29 {+/-} 8.9 months. There were significant reductions in total grey matter and white matter volumes, together with volume loss in the bilateral hippocampus, bilateral parahippocampal gyrus and right entorhinal cortex (pFDR[≤]0.05) over the follow up period. Volume loss in the right parahippocampal gyrus was associated with increase in aggression (pFDR<0.01), and CBF increase in the right hippocampus and parahippocampal gyrus was associated with worsening in agitation (pFDR<0.01). Furthermore, lower baseline volume of the right entorhinal cortex was associated with increase in aggression, and lower baseline CBF in the right hippocampus and parahippocampal gyrus was associated with worsening of agitation (pFDR<0.05). Our findings highlight the importance of identifying NPS in individuals with MCI and the role of structural and perfusion MRI measures as promising biomarkers for monitoring disease progression and improving risk prediction of conversion to dementia.
Ling, Y.; Tsung, S.; Alarid-Escudero, F.; Goldhaber-Fiebert, J. D.; Brandeau, M. L.; Reitsma, M. B.
Show abstract
People experiencing homelessness (PEH) have substantially lower quality-adjusted life expectancy and substantially higher healthcare spending compared to housed populations. California has the largest population of PEH in the United States, and an estimated ten percent of PEH have co-occurring opioid use disorder (OUD). We developed a lifetime, individual-level microsimulation of the 18,928 PEH with OUD in California to evaluate the health and economic effects of expanding access to housing through a Housing First approach, which does not require OUD treatment or abstinence. Increases in permanent housing and shelter access were based on Los Angeles County's Measure A goals, and we accounted for both housing and healthcare costs. Expanding housing access increased mean age at death by 0.34 years, reduced fatal overdoses for the cohort by 39, and produced an incremental cost-effectiveness ratio of $31,083 per quality-adjusted life-year gained, well below conventional willingness-to-pay thresholds. We provide evidence that housing is a high-value health intervention for PEH with OUD in California. Policymakers should support Housing First programs as part of integrated responses to the intertwined crises of homelessness and opioid use.
Kagawa, E.; Kato, M.; Hidaka, T.; Kunita, E.; Hironobe, N.; Matsui, S.; Kobayashi, Y.; Uotani, Y.; Teramoto, T.; Yoshii, K.; Okumoto, T.; Tani, M.; Dote, K.
Show abstract
Background: Artificial intelligence (AI) may improve electrocardiographic differentiation of Takotsubo syndrome from acute anterior ST-segment elevation myocardial infarction (STEMI), but its impact on human diagnostic performance, the added value of visual explanations, and differences by reader background and experience remain uncertain. Methods: Using 12-lead electrocardiogram (ECG) images from two hospitals, we developed and externally validated an ensemble of three deep learning models after self-supervised pretraining and fine-tuning. The diagnostic task used ECGs recorded within 14 days after symptom onset, regardless of reperfusion status. We conducted a four-phase reader study. Phases 1-3 compared unaided interpretation, a reference sheet, and additional AI assistance. In Phase 4, readers evaluated a new ECG set after randomization to reference-sheet support alone, with AI predictions, or with AI predictions plus gradient-weighted class activation mapping (Grad-CAM) heatmaps. Diagnostic correctness was analyzed using generalized linear mixed-effects models. Exploratory analyses examined differences by cardiologist status and experience. Results: Development and external-validation cohorts, including healthy controls, comprised 832 and 95 participants, respectively. For Takotsubo-STEMI discrimination, external-validation accuracy was 0.914, with areas under the receiver operating characteristic and precision-recall curves of 0.970 and 0.984. Seventy healthcare providers completed the reader study. The estimated probability of correct diagnosis increased from 0.632 without assistance to 0.864 with AI assistance (odds ratio [OR], 3.71; P<0.001). In randomized Phase 4, AI assistance improved accuracy vs reference-sheet support alone (OR, 4.28; 95% CI, 3.08-5.95; P<0.001), and Grad-CAM provided additional benefit (OR, 1.53; 95% CI, 1.04-2.25; P=0.033). Exploratory analyses showed greater improvement among less experienced readers; with AI assistance, non-cardiologists outperformed cardiologists. AI-assisted readers remained less accurate overall than the model alone. Conclusions: An ECG image-based ensemble differentiated Takotsubo syndrome from anterior STEMI and improved reader accuracy. Visual explanations provided additional benefit, while experience-related differences highlight the importance of effective human-AI interaction.
Zhigulev, A.; Poujade, F.-A.; Gryzunov, N.; Buyan, A.; Li, H.; Bergman, O.; Liu, M.; Pradhananga, S.; Franco-Cereceda, A.; Swoboda, P.; Kulakovskiy, I. V.; Eriksson, P.; Sahlen, P.; Björck, H. M.
Show abstract
BACKGROUND: The chromosome 9p21.3 locus is among the strongest genetic risk loci for cardiovascular disease (CVD), yet the regulatory mechanisms underlying this association remain poorly understood. Previous vascular disease studies focused predominantly on the role of this locus in smooth muscle cell (SMC) dysfunction, whereas its potential role in endothelial cells (ECs) remains less well characterized. Here, we sought to identify candidate regulatory variants at 9p21.3, their distal target genes, and effector transcripts that contribute to CVD-related phenotypes in ECs. METHODS: Promoter capture Hi-C was used to map chromatin interactions at 9p21.3 in primary human ascending aortic SMCs (n=6), supplementing the previously obtained ECs (n=16) dataset. To support the regulatory relevance of these interactions, expression quantitative trait locus analysis was performed in 136 vascular biopsies. Variant effects were prioritized in silico using AlphaGenome. The regulatory activity of the prioritized region was validated using CRISPR activation. RESULTS: Chromatin interaction mapping identified a complex regulatory architecture at the 9p21.3 locus involving 19 putative regulatory elements and 15 genes across ECs and SMCs. MIR31HG showed the most recurrent distal promoter connectivity with the 9p21.3 region. Consistent with the chromatin interaction map, expression of MIR31HG transcripts was associated with genetic variation at the 9p21.3 locus. In ECs, the lncRNA MIR31HG, but not miR-31, was responsive to hemodynamic flow. Depletion of MIR31HG in ECs was associated with endothelial activation. Epigenetic tracks and deep learning predictions highlighted an active enhancer element at 9p21.3, which was subsequently validated by CRISPR activation and increased MIR31HG expression, confirming its capacity to regulate MIR31HG over a large genomic distance of approximately 538 kb. CONCLUSIONS: Altogether, these data implicate lncRNA MIR31HG as a likely effector of the 9p21.3 CVD locus, linking the cis-regulatory element and heritable genetic variation in the respective region to endothelial activation through long-range enhancer regulation.
Bou Ali, M. N.; Lassoued, N.; Ben Amor, B.; Nasr, S.; Aouididi, M. A.
Show abstract
Background. In people with type 2 diabetes (T2D), postprandial glucose weighs on glycaemic control, and carbohydrate counting is difficult. Patients can describe their meal but rarely quantify it. Objective. To assess (1) whether a small open large language model (LLM), fine-tuned on public data, estimates the carbohydrate content of a meal described in words, and (2) what this carbohydrate adds to the prediction of the postprandial glucose rise in two public cohorts. Methods. Step 1: the Qwen3-4B LLM was fine-tuned with QLoRA on Nutrition5k dishes and compared with an automatically matched USDA food composition table and with Qwen3-4B and Qwen3-8B used without training, on test dishes (internal control) and on 100 real meals of ShanghaiT2DM patients, whose reference carbohydrate was computed from the recorded quantities. Step 2: the maximum glucose rise between 0 and 120 minutes ({Delta}Gmax) was predicted by personalised gradient boosting (M2), compared with a spline linear regression and with the Chronos-2 foundation model, using only information available at meal start and person-level cross-validation, in ShanghaiT2DM (people with T2D recruited in two Shanghai hospitals; carbohydrate estimated by the LLM) and CGMacros (free-living adults in Santa Barbara, California; carbohydrate entered in the dataset, no text). Step 3: the model was transferred between datasets. 95% confidence intervals (CI) were obtained by cluster BCa bootstrap (2,000 resamples). Results. On the 100 real meals, the carbohydrate mean absolute error (MAE) of the fine-tuned model was 12.2 [9.9, 15.2] g, almost half that of the table (22.6 g; paired difference -10.5 [-16.8, -6.0] g) and nearly a third of that of the same model without fine-tuning (34.3 g), with a bias of -10.3 g; on test dishes it was 2.2 [2.0, 2.5] g. Personalised gradient boosting (M2) reached a correlation between predicted and observed rises of r = 0.63 [0.55, 0.69] in CGMacros (1,286 meals, 45 adults) and 0.49 [0.42, 0.54] in ShanghaiT2DM (2,339 meals, 99 people with T2D), about twice the linear reference (r = 0.32 and 0.23); for a rise above 50 mg/dL, its area under the ROC curve (AUROC) was 0.80 and 0.74, versus 0.64 and 0.68 with pre-meal glucose alone. Adding carbohydrate to pre-meal glucose changed r by 0.095 [0.068, 0.128] in CGMacros and by 0.004 [-0.010, 0.017] in ShanghaiT2DM; difference between datasets: 0.091 [0.060, 0.126]. Applied from one dataset to the other, the model lost correlation ({Delta}r = -0.11 [-0.17, -0.05] towards ShanghaiT2DM and -0.12 [-0.19, -0.04] towards CGMacros). The full "text [->] LLM [->] glucose" chain was evaluated in ShanghaiT2DM only. Conclusion. An open 4-billion-parameter LLM, fine-tuned with QLoRA on public data and usable offline, almost halves the error of carbohydrate counting for real meals described in words compared with an automatically matched table, with an underestimation of large meals. Personalised gradient boosting, based on pre-meal glucose and the person's usual response, doubles the correlation of a linear reference for predicting the glucose rise in both cohorts. Carbohydrate improved this prediction in CGMacros but not in ShanghaiT2DM; the two datasets differ in many respects, and the cause of this gap is not established.
Arojo, O.; Jones, A. M.
Show abstract
Background: Epidermal Growth Factor Receptor (EGFR) tyrosine kinase inhibitors (TKIs) are widely used for EGFR-mutated non-small cell lung cancer (NSCLC) but have distinct toxicity profiles across generations. Aims: To investigate relationships between EGFR-TKI physicochemical, pharmacokinetic and pharmacological properties and drug-specific adverse drug reactions (ADRs) using national and global pharmacovigilance data. Methods: ADRs for erlotinib, afatinib and osimertinib were extracted from the UK MHRA Yellow Card scheme and WHO VigiAccess. Disproportionality analysis used reporting odds ratios (RORs). Signals were defined as ROR >2, lower 95% confidence interval >1 and >3 cases. Drug properties were curated from ChEMBL, emc, SmPC and EMA sources. Results: VigiAccess contained 88,662 individual case safety reports (163,619 ADRs) and Yellow Card contained 752 reports (1,814 ADRs). Four System Organ Class (SOC) signals were identified in VigiAccess and three in Yellow Card, with five additional Yellow Card Preferred Term signals. Eye disorders showed the strongest consistent association with erlotinib (VigiAccess ROR 4.33; Yellow Card ROR 4.18). Osimertinib was associated with pneumonitis (ROR 7.67) and pulmonary embolism (ROR 3.94). Afatinib showed psychiatric disorder signals in both databases (ROR 3.2-3.3) and diarrhoea (ROR 2.79). Neoplasm-related signals likely reflected confounding by indication. Conclusions: EGFR-TKIs demonstrate generation-specific pharmacovigilance profiles. Integration with drug properties supports mechanistic plausibility for several signals, while the atypical afatinib psychiatric association warrants further investigation. Findings are hypothesis-generating and require confirmation in controlled observational or mechanistic studies.